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I. INTRODUCTION 



The Dirac-Kepler problem in D + 1 dimensions has been treated in several ways: by power 
series^, radial SUSY QM^, supersymmetry generated by the Lippmann- Johnson operator- 
and an su{l, 1) approach^. On the other hand, the Dirac radial equations with vector and 
scalar Coulomb-type potentials were studied by series power-, SUSY QM^ and intertwining 
operators^. 

Joseph was the first who studied the Dirac equation with vector potential in D spatial 
dimensions^. Also, the Lippmann- Johnson operator was introduced to generate the super- 
symmetric charges^. With Coulomb-type scalar and vector potentials, the problem was 
solved in analytical way by reducing the radial Dirac equations to the differential equations 
satisfied by confluent hypergeometric functionsi^i^. The energy spectrum and the SUSY 
ground state of this problem were found by an sm(1, 1) algebraic approach^. However, to 
our knowledge the Lippmann- Johnson operator in general dimensions with Coulomb-type 
scalar and vector potentials has not been reported. Therefore, the supersymmetry generated 
from this constant of motion remains untreated. 

The purpose of this paper is to construct the supersymmetry charges for the Dirac-Kepler 
problem plus a Coulomb scalar potential from the generalized Lippmann- Johnson operator 
in D dimensions. The fact that one of the supercharges annihilates the SUSY ground state 
leads us to find the energy spectrum for discrete excited states and to obtain the radial 
differential equations for the SUSY ground state. By performing a similarity transformation 
to the radial Lippmann- Johnson operator we find the SUSY ground state. Also, we show 
that the radial part of the Lippmann- Johnson operator is reduced to that reported in Ref. 
[7 for the three-dimensional space and obtained by intertwining considerations. 

II. THE RELATIVISTIC DIRAC EQUATION IN D + 1 DIMENSIONS AND 
THE GENERALIZED LIPPMANN-JOHNSON OPERATOR 

The Dirac equation in D + 1 dimensions for a central field can be written with {h = c = 1) 

asr- 

Hm = |7Vy + 7° + K (r)) + K (r) }^ = J = 1, 2, D (1) 
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where summation over repeated index is assumed, m is the mass of the particle, Vg and 
Vv are the spherically symmetric scalar and vector potentials, respectively, and the Dirac 
matrices in D dimensions, 7-', satisfy {7^,7*^} = 277^^^ with 

{Sf"" if u = 0, 
(2) 
-5^^ if /X ^ 0. 

In D spatial dimensions, the orbital angular momentum operator Lab and the total an- 
gular momentum Jab are defined as 

Lab = ixadb - ixbda (3) 

and 

Jab = Lab + ^7^7', (4) 

respectively. In the case of spherically symmetric potentials, the total angular momentum 
operator and the spin-orbit operator 

KD = -A\Y.^^l'Lab + \{D-l)\ (5) 

commute with the Dirac Hamiltonian. For a given total angular momentum j, the eigenval- 
ues of Kd are = ± (j + 1/2), where the minus sign is for aligned spin j = ^ + |, and the 
plus sign is for unaligned spin £ — i^. 

We find that if the scalar and vector potentials are given by V^^ = ^ and = ^ then, 
the matrix hermitian operator 

B = -iKdi''^^ [H - 7°m) + 7''+S°f ^ («.m + a^H) (6) 

is a constant of motion (see appendix A), where the pseudoscalar 7^+^^ is reduced to the 
matrix 7^ in (3 + 1) dimensions and satisfies (7"^^^)^ = 7^^"*^, (7^^"*^)^ = 1 and {7"^^^, 7^} = 
0. Also, it can be shown that this operator anticommutes with the Dirac operator Kd. 
Therefore, the operator B satisfies 

B^!,^ = -6v[/_,^, (7) 
where h is an undetermined constant. In fact, B is the generalization of the Lippmann- 



Johnson operator-i^ and is reduced to that given in Ref. 



7 for D = 3, and to that for D 



dimensions in absence of the scalar potential Vg reported in Ref. 



In this way, for odd or even dimensions, we write the eigenstates of equation ([T]) as 

^^D=r ^ \ ... e (8 



being G'1^2('") the radial functions, and (^d) the hyperspherical harmonic 

functions coupled with the angular momentum j^. We consider that 

l^^Vf^^.n = -^-^n, (9) 

which is the generalization of the three-dimensional equation {a ■ r) \E'k = — ^-k-*^- By 
defining the spinors 

= I "^^"^^"^ I e-^^*, = I I e-^^* (10) 

and from the results of Appendix B, the explicit form of the operator B acting on a general 
eigenstate of the Dirac Hamiltonian is 

+ a,7°(m + K + 7°K)$«,}. (11) 
III. SUSY QM AND THE ENERGY SPECTRUM 



Based on Refs. l3lJ9lJl3l. we define the supersymmetric generator 



.00. 

Q=\ , 12 

BO 



which satisfies {Q, Q} = and 



n^{Q,Q^}=\^ °|, (13) 
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with "H the supersymmetric Hamiltonian. 

In order to obtain the energy spectrum for the Dirac Hamiltonian H, we consider the 
results given in Appendix C, from which 

B^ = {a^m + asHf + {H^ - m^) . (14) 
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By considering equation ([7]), we have the eigenvalue equation = ^^^k^,- Thus 

^2 = ^a^m + asEf + (E^ - m^) . (15) 
Since the supersymmetric ground state, "^susy must satisfy the condition 

n^lusY = 0, (16) 

it follows that the ground state energy eigenvalue, Eq, is obtained from equation f[T^ by 
setting 6 = 0. In this way 
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E.-'ni-^±Jl4^] -l^]y (17) 



where 7^ = /tfj + — a^. For the excited states of the Hamiltonian H, we perform the 
change 7 — )■ 7 + n, where n = 0,1, 2, 3... is the radial quantum number. Then 



K_ a.a. ±\l(-r-T^X-(4^^P^]^ (18) 



m Q;2 + (7 + n)2 V V'^? + (t + "')^/ + (7 + ^)^/ ' 

which is in accordance to that obtained from an analytical^*^ or su{l, 1) algebraic approach^. 

The eigenstates \E'k^ and "^-kd transformed into each other by the operator B (equa- 
tion ([7])) and both are eigenfuntions of the operator with the same eigenvalue. Therefore, 
the supersymmetric eigenstates can be written as 

'isusY = I I (19) 

Considering equations ( [T3l) and ( IT6l) . the components of the sypersymmetric ground state 
must satisfy 

B^o±.o = 0- (20) 

In order to solve this equation for +kd (the solution for the other sign can be obtained by 
equation ([7])) we consider the expression (lB3p . from which 



>■ ' ' -•G'S.X- 



( I y (21) 



Thus, the radial components of the SUSY ground state {Gq I^ and Gq I^) satisfy 

(22) 




where 




K]j dr r 



(23) 



dr r KD 



e± = K£, + a^a^l and a± = (a^ ± a^). The matrix operator L^, obtained by means of 
SUSY QM, is reduced to the three-dimensional operator L reported in Ref. 13, which has 
been constructed by imposing an intertwining relation between the corresponding radial 
Dirac Hamiltonian and L. 

In order to find the expression for the radial components of the SUSY ground state, we 
define 




^(2) 

Therefore, equation fl23|) is rewritten as 

d , 1 / a+ I W GqId \ '^i^D 



(24) 




which can be easily solved by diagonalizing the matrix of the factor 1/r. For this purpose, 
we perform the transformation 




Thus, from equation fl25l) . we obtain 

'^i) + '^s I otvi^D — Olsl 

The unnormalized solutions for these differential equations are given by 



(26) 




)l 






1^(2) 1 



(27) 



Poln ^ exp ( 2——^ i^^'^D + a.7) ) , (28) 

fI^L = ( 2 ^ „,2 («^^^ ~ ««7) ) • (29) 
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Since -Pq diverges at r = 0, it is not a physically acceptable solution. Hence, the radial 
spinor for the supersymmetric ground state is 



^susY - I ^ / ° X I • (30) 



For the case = 0, ipsusY ^ normalizable solution only for < which is in accordance 
to the results reported in Ref. y for 3 + 1 dimensions. 

Having obtained the SUSY ground state, the explicit form of the eigenf unctions cor- 
responding to higher SUSY energy levels should be determined by solving the eigenvalue 
equation 'H'^susy = b'^'^susY which is equivalent to find the solutions of the equation 
S^^'k^ = fe^^E'reo- Nevertheless, this problem is much more complicated to solve than the 
original Dirac eigenvalue equation. 



IV. CONCLUDING REMARKS 

We treated the Dirac-Kepler problem plus Coulomb-type scalar potential by generalizing 
the Lippmann- Johnson operator in general dimensions. This operator allowed us to con- 
struct the supersymmetric charges from which we found the ground state energy spectrum 
and the spectrum energy for discrete excited states. The action of the Lippmann- Johnson 
operator on the ground state of the Dirac equation leads to the radial operator Ljy, equation 
fl23|) . which generalizes to D-dimensions that reported in Ref. Q for the three-dimensional 
space. The ground state obtained in this paper is in full agreement with those reported in 



Ref. 



With the generalized 



141 and it reduces to that in three dimensions reported in Ref. 
Lippmann- Johnson operator reported in the present paper, we can construct an 5*0 (4) sym- 
metry treatement, similar to that given for the Kepler- Coulomb problem in three spatial 
dimensional^, which is a work in progress. 



Appendix A: Calculation of [B,H] = 

For an arbitrary radial function /(r), we show that [Ko,f{r)] = 0. Moreover, we find 
that 



Hence, 



-27^+S° (m + K) 



X" 



2i 



(Al) 
(A2) 

(A3) 



[5, H] = -iKd [H, 7^+1] {H - 7°m) 
+ «mir^7^+i [iJ,7°] 



+ 



i7,7^+SV 



: -2^ir^7^+S°{(m + V^s) - 7°^) 
0. 



(A4) 



Appendix B: Calculation of -B^I/k^ 

Considering the properties of the operator Kd, the algebra satisfied by the matrices 7* 
and from equations ([9]) and (fTOj) . we obtain 



b i f ot^ D-1 
r r \ 2 



+ (K-7V,) e.,} 



(Bl) 
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and 



+ [q;,7° {m + Vs + 7°K) + ma^] } . (B2) 

Thus, the exphcit form of the Lippmann- Johnson operator B acting on an eigenstate 
of the Hamiltonian H is 

= -r-^l {a^dr - kdK) Oku + {^^D (^r + ^7°) + ma.,^ 

+ a,7°(m + K + 7°K)$,,}. (B3) 



Appendix C: Calculation of 

With the definitions A\ — H — 7°m and A2 = a^m + a^if, we find the following commu- 
tation relations 



Since 



[^1,^2] 
[Ai,7^^^] 
[^2,7^-^^] 



A,7V- 
r 



— 2mQ;57V, 
-27^+S°K, 

2^7° 



^2, 7 V 



-i^D + 2Ky-, 



2iK7°^i^ + 2a. (m + 7 



a; 



-f^D, 7°7*— = and {i^D, 7'°'^^} = 0, we finally obtain 



(CI) 
(C2) 
(C3) 

(C4) 
(C5) 



52 = B^B = + - ai7°y-A2XD 

r 



a; 



+ a27V- Ai^D 

r 



= + + 2m (K + l^'Vs) 7°i^^ 



(C6) 
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